Counting statistics for electron capture in a dynamic quantum dot 
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We report non-invasive single-charge detection of the full probability distribution P n of the ini- 
tialization of a quantum dot with n electrons for rapid decoupling from an electron reservoir. We 
analyze the data in the context of a model for sequential tunneling pinch-off, which has generic 
solutions corresponding to two opposing mechanisms. One limit considers sequential "freeze out" of 
an adiabatically evolving grand canonical distribution, the other one is an athermal limit equivalent 
to the solution of a generalized decay cascade model. We identify the athermal capturing mecha- 
nism in our sample, testifying to the high precision of our combined theoretical and experimental 
methods. The distinction between the capturing mechanisms allows to derive efficient experimental 
strategies for improving the initialization. 

PACS numbers: 



The fast formation of quantum dots (QDs) out of a 
two-dimensional electron system (2DES) constitutes an 
open problem within the field of nanoscale electronics [l[ . 
The initialization process in these dynamic QDs is a key 
ingredient in, e.g., devices for quantum information pro- 
cessing 0, S] , single-electron current sources , or na- 
noelectronic circuits 0, S] ■ The outcome of the initial- 
ization is characterized by a probability distribution P n 
for trapping exactly n electrons in the QD. The goal is 
to attain a predictable low dispersion distribution thus 
making dynamic QDs reliable and reproducible sources 
of electrons on demand. Deviations from this ideal case 
may be caused, for instan ce, by backtunneling 0-11 or 



non-adiabatic excitations 12-141. 



A decay cascade model |rij has been proposed recently 
to predict P n in dynamic QDs. It has become popular 
for benchmarking QD-based current sources [15l4l8| by 
extracting P n from the average current as function of 
QD energy using the model. Alternative mechanisms, 
such as sudden decoupling from thermal equilibrium [l9| 
have been proposed. Experimental distinction between 
the capturing mechanisms is the key towards systematic 
improvement of the initialization precision. So far P n 
has not been measured with sufficient accuracy to allow 
this distinction. The first two cumulants of P n have been 
extracted from current and noise measurements [20|-|22j. 
Furthermore, single charge detection [l9|, 23 1 has been 
used to determine partial information on the distribution 
of QD population-depopulation events. 

Here we present non-invasive charge detection to mea- 
sure the full probability distribution P n of capturing 
n electrons in a dynamic QD. Considering the integer 
charge on the QD to be the only degree of freedom out of 
equilibrium, we derive theoretically two generic limits for 
P n : a (generalized) frozen grand canonical distribution 




Figure 1: (Color online) (a) False-color electron microscopy 
image of a device. The upper half shows the semiconducting 
part consisting of an 800 nm wide channel (light blue) crossed 
by topgates (yellow) . The QD is formed by the leftmost group 
of topgates, between source (S) and drain (D). The light-gray 
parts in the lower part form the SETs labeled Detl (red) and 
Det2 (blue), respectively, (b) Detector signals of the charge 
transfer sequence, as explained in the main text. 



and a rate-driven athermal limit (generalizing the decay 
cascade model ll|). In the generalized form both limits 
may be hard to distinguish experimentally for P n ~ 1, 
which is the relevant regime for most applications. Yet, 
our experimental data for P n allow to distinguish both 
limits and to conclude that the dynamic QD initializa- 
tion is consistent with the athermal distribution. Based 
on these findings, strategies for optimum high fidelity ini- 
tialization are presented. 

The device under investigation is shown in Fig. [TJa) . 
The top half of the false-colored image shows four QD 
structures in series consisting each of three gates cross- 
ing a 2DES within a AlGaAs/GaAs heterostructure. The 



2DES is located 90 nm below the surface, the wet- 
etched channel is 800 nm wide. Similar QD struc- 
tures have previously been used as single-electron cur- 
rent sources [24| . We use the left QD as the dynamic 
QD, which captures electrons from source (S) and after- 
wards emits them to the node (dotted region) for charge 
detection. The voltage on one gate of the rightmost QD, 
Vbarricr, controls the transparency of the node to the drain 
lead (D). All other gates are grounded and do not affect 
the circuit. Close to the node two single-electron tran- 
sistors (SET) based on Al-AlO^-Al tunnel junctions are 
placed as charge detectors (Detl/Det2). To increase the 
coupling between the potential of the node and the metal- 
lic detectors both are capacitively coupled by H-shaped 
floating gates (FG). Correlating the detector signals al- 
lows to distinguish the electron signal from background 
charge fluctuations, as both detectors are coupled to the 
same island. All measurements are performed in a dry di- 
lution cryostat at nominal temperature of about 25 mK. 

Fig. [H[a) shows the sequence of QD initialization and 
charge transfer to the node schematically. An isolated 
QD (ii) is formed between the two left-most gates by 
applying sufficiently negative dc voltages, Vgs and Vgd 
(see Fig. [TJ. Initialization (i) of the QD is achieved by 
applying the first half cycle of a sinusoidal pulse super- 
imposed onto the source gate, so that the source barrier 
becomes transparent. During the subsequent rise of the 
source barrier a certain charge state of the isolated QD 
(ii) with n electrons is established with probability P n . 
In the second half cycle of the sinusoidal pulse (iii) the 
left barrier is raised further. As the source barrier also 
couples to the QD potential, i.e. acts as plunger gate [2~ij |. 
one can ensure complete unloading of the QD charge to 
the node where it can be detected non-invasively. 

The cycle is repeated three times during which charges 
accumulate on the node. Afterwards, opening the right 
barrier resets the node charge, before the next three cy- 
cles start. Example traces of the SET signals are shown 
in Fig. [lib). The bold black vertical lines represent resets 
of the node's charge state. The stochastic nature of this 
process is represented by the different initial states after 
each reset. Each thin green vertical line represents a com- 
bined charge capture and transfer pulse. Due to the dis- 
creteness of node charges, we can assign levels (horizontal 
lines) derived from a histogram of the detector trace to 
each interval between pulses. These are then compared 
to extract the number of captured electrons in this spe- 
cific cycle (see Appendix) . Due to the limited bandwidth 
of the SET detectors (/b w ~ 600 Hz) the pulses are de- 
layed by 40 ms each. The pulse itself consists of a single 
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period of a sine with frequency / pu 



40 MHz. 



The voltage Vgd allows to adjust the depth of the 
QD potential and thereby to tune the average number 
(n) = J2n of captured electrons [25j- Fig.[3]shows P n 
as measured by charge detection as function of Vgd • The 
probabilities of charging the QD with up to 4 electrons 
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Figure 2: (Color online) (a) Sequence of schematic potential 
landscapes during capturing from source (i), isolation (ii) and 
emission to drain (iii), respectively. Shading indicates tunnel- 
ing rates F n , related to area under curve, (b) Evolution of 
the tunneling rate T n and energy p n of a particular charge 
state for two different Vgd voltage settings. 



are well resolved. For n — 1 the initialization accuracy 
reaches 99.1 % for V" G d ~ -192.5 mV. The probability 
to charge the QD with 4 electrons with one initializa- 
tion pulse approaches 80 % for Vgd ~ — 168 mV. We 
will analyse the V"gd dependence of P n further below. In 
the following the theoretical framework is established al- 
lowing later to relate the measured distribution to the 
underlying capturing mechanisms. 

Treating tunneling perturbatively in a Markov ap- 
proximation 0, [53] , the "disequilibration" and eventual 
freezing of P n {t) can be described by a general master 
equation: 



Pn (t) = P„_J (t)W+_ 1 (t) - P n (t)W~ (i)+ 
Pn+l{t)W- +1 {t)-P n {t)W+{t), 



(1) 



where are the instantaneous rates for adding (+) or 
removing (— ) an electron to/from the QD, averaged over 
all degrees of freedom except n. The balance between 
adding and removing defines (possibly non-equilibrium) 
electrochemical potentials /z„(t) of a QD state with n 
electrons: 



W-(f)/Wtf_i(i). 



(2) 



with f3 = kT being the product of temperature and 
Boltzmann's constant. If the time-dependence of rates 
W^(t) is quasistatic, then Eq. ([2]) is the expression of 
thermodynamic detailed balance and fj, n = Q n — 
will be set by the differences of the thermodynamic po- 
tentials 27] Q n = J- n + E n — nEp, where J- n is the 



3 



canonical free energy of the internal degrees of freedom 
on the QD, E n is the electrostatic interaction energy and 
Ep = is the Fermi level in the source lead. Defining 
the total rate for charge exchange in no?i — 1 transition 
as r„(i) = W~(t) + W+_ x (i) allows us to write Eq. {J) 
in the form 

Pn{t) = -r„ [/(^ n )P„(t) - /(/X„)P n _l(t)] + 

r n+1 [/( Mn+1 )p„ +1 (t) - /( Mn+1 )p n (t)] , (3) 

where f(x) = 1/(1 + e^ x ) is the Fermi distribution in the 
source and /=! — /. 
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Figure 3: (color online) Probabilities to capture n = 0...4 
electrons per cycle derived by counting as a function of Vgd. 
The error bars indicate the 95 % confidence interval. 

We will now apply Eq. Q to the decoupling process 
sketched in Fig. [5Ja-i)-(a-ii). The end of the decoupling 
stage (ii) is characterized by T n (tf) = 0. The exact time- 
dependencies of {fJ- n (t)} and {T n (t)} between an initial 
time moment to and t / are impossible to predict without 
a more specific microscopic model for the QD. However, 
two sets of critical parameters can be unambiguously de- 
fined. For each transition between n — 1 and n charges on 
the dot we identify two time moments: one for the onset 
of backtunneling, i„, and the other one for decoupling 
(detailed balance breakdown), t°, described in Fig.[2](b). 
It shows schematically the evolution of T n (blue) and 
fi n (green) with time for the initialization process with 
two different values of Vgd • The time for onset of back- 
tunneling is defined by the crossing of the Fermi level, 
Hn{t b n ) = 0, while the decoupling time is set by the av- 
erage number of remaining tunneling events being of or- 
der one, yj T n (t)dt = 1 (shaded area under the curve). 
Positive charging energies and the raising bottom of the 
confining potential well imply t b n+1 < t„ while in the 
relevant regime we expect the higher n states to be less 
stable, T n+ i(t) > T n (t) (at least for t > t b n ) and therefore 
generally t c n+1 > t°. Sufficiently long equilibration before 
the decoupling implies that t c n is well-defined for all n for 
which the Coulomb blockade holds. Fluctuations in the 
capture probability (i.e. P n ) will be strongest for n that 
have t£ and t h n close to each other [28 1. 



We proceed by solving Eq. (J3]) for P„ = P n (t tf) 
in two limits (thermal versus athermal) that correspond 
to opposing physical mechanisms of charge capture. We 
assume that initially the charge on the QD is equili- 
brated, corresponding to a grand canonical distribution, 
P n (to) e ~P n -n.{t ) _ As long as T n {t) remain sufficiently 
large, the solution P n (t) at t > t closely follows the in- 
stantaneous equilibrium. This can be seen directly from 
Eq. (J3j) : large T n pin the terms in square brackets to zero 
and the evolving distribution of n obeys detailed balance 
adiabatically, 



P n (t)« e ft-+iWp n+1 (t). 



(4) 



In deriving the thermal limit we consider sudden decou- 
pling, i.e. r n+ i(i) dropping to so fast that Eq. (|4]) 
must hold up to t = t„ +1 but once t > t° +1 the r.h.s. 
of Eq. ([3]) is effectively zero and P n {t) freezes (i.e. re- 
mains constant). With this sudden approximation, the 
asymptotic value P„ = Pn(*n+i) i s se * by a "curtailed" 
grand canonical distribution that excludes already frozen 
charge states with n' < n but is normalized over the re- 
maining states with n' > n that keep being connected 
until t c n+2 . This gives P n = (1 - Em=o P ™) Z n+i or ex- 
plicitly 



(5) 



where Z n = l+E^nUtn^ 1 ^ = l+e-^«)(l+ 
e -/3M>i+i(*„)(T -)- _ )) includes the electrochemical po- 
tentials jj, n i of states n' > n taken at the decoupling 
moment t c n of the state n. The assumption of well- 
pronounced Coulomb blockade implies that the addition 
energy, A/i(t) = fx„(t) — fx n -i(t) remains large compared 
to temperature. Thus there exists sufficently large St 
such that (3Afj,(t) > 5p for all relevant n and t. For 
<5t > 1, a further approximation, Z n w 1 + e - ^"'*™), 
results in at most e~ &T relative error for each P„, leading 
to a simple expression for the low-dispersion limit of the 
sudden decoupling mechanism, 



^«=/[/W*n+l)] II f^rn(t C r , 



(6) 



The distribution ([6]) is determined by a set of dimension- 
less numbers jl n = j3/j, n (t^ l ) and is narrowly dispersed if 
• • • > /Wi > Mm > ■ ■ ■■ 

Now we consider the athermal limit. At sufficiently 
low temperatures the timescale for f(p n (t)) switching 
between loading (w 1) and unloading (w 0) may be- 
come much shorter than the time-scale of reducing r„ (t) . 
Assessing this gradual decoupling limit amounts to re- 
placing the Fermi functions in Eq.([3|) by sharp steps, 
f{Hn) — > Q(t — t b n ). Starting with a sharp initial equi- 
librium free of thermal fluctuations, P n {to) — <5 n .Ar with 



4 



. , the system of equations Eq.Q is re- 
duced to a set of decay cascade equations: 



C +1 < <o < 



P n (t) = 




t < t 



n+l) 
-1 < t < t b n , 



t" < t 



(7) 



with -Po(^) set by normalization. These equations gener 



alize the decay cascade model [VL| to distinct f„'s. One 



can show that a universal solution to Eqs. (0 indepen- 
dent of the specific shape of r„(t)'s time-dependence (llj . 



TV 



Pn 



-x, 



n a 

j=n+l 



remains valid in the limit of . 



> X i+1 >Xj » 



(8) 



with 



appropriately generalized integrated decay rates, X n = 
Jt{T n (t)dt T ±° fl f W-(t)dt- 

Next we compare the two theoretical limits (j6|) and (jHJ) 
with our measurement of P„, as shown in Fig. |4] In the 
main diagram the experimental data for Pi are plotted 
as symbols as function of Vgd ■ As indicated in Fig. [2][b) , 
Vgd controls £t n and X n . We assume fx n to be linearly 
dependent on Vgd while typically decay rates depend 
exponentially on the gate voltage. Thus we use 



lnX, 



(9a) 
(9b) 



incorporating Vqd dependence in Eqs. ([5]) and (J5J using 
a^n and as fit parameters for the thermal model 

and ax, n and Ax, n for the athermal model. The red 
and black solid lines in Fig. [4] show the result of fit- 
ting the thermal [Eq.(j6|) and (JHK)] and athermal distribu- 
tion [Eq.([SJ and (J5b)] to the experimental data around 
Pi » 1. The blue, black and green data points corre- 
spond to Pq, Pi and P2 respectively, shown on a loga- 
rithmic scale in the inset. Only OhI and 1^2 transi- 
tions have been considered since P„ with n > 2 do not 
contribute for Vgd < — 186mV (see Fig. |3J). The er- 
ror bars on the data indicate 95 % confidence intervals 
for estimating probability from the binomial statistics of 
direct counting [30] (see Appendix). The thermal limit 
(red solid line) deviates beyond confidence interval for 
voltages Vgd > — 190 mV for Pi on linear scale (see red 
arrow). Note also that the uncertainties for the model 
parameters turn out larger for the thermal model (see 
Table fl} . Even stronger deviations can be observed for 
P2 on the logarithmic scale. Here the thermal model 
predicts a linear characteristic, while the apparent non- 
linearity is reproduced well by the athermal limit. In fact, 
the non-linear characteristic of P2 is a direct signature 
of simultaneous variation in tunneling rates and energy 



Table I: Best fit parameters. 





a M /l,l/mV 1 


otf,/x,2/mV 1 


A M /x,i 


A tl/X ,2 


Thermal 


-0.293 ± 0.003 


-0.983 ± 0.023 


-60.8 ± 0.6 


-182.1 ± 4.2 


Athermal 


0.261 ± 0.003 


0.385 ± 0.009 


-54.5 ± 0.6 


-71.5 ± 1.6 
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Figure 4: (color online). Measured Pi (points) as function 
of Vgd, compared to theoretical fits to thermal and athermal 
limit (red and black solid lines, respectively). The error bars 
indicate 95 % confidence intervals. The blue, black and green 
data points correspond to Po, Pi and P2, respectively, shown 
on a logarithmic scale in the inset. 



during the capture process. Thus we conclude that our 
measured distribution is consistent with the generalized 
decay cascade model in the low- noise limit. 

The generalized decay cascade limit implies that low- 
ering the lead temperature will not increase the initial- 
ization precision, which has indeed been found in sur- 
face acoustic wave driven devices [3l[ . This feature alone 
would, however, not suffice to exclude the thermal distri- 
bution as in the past this saturation has been related to 
rf-heating induced by the modulation 31[. Our results 
indicate a path for further improvements for dynamic QD 
initialization: it may be achieved by increasing the sepa- 
ration of decay steps (X n /X n _i) either by a large decay 
rate ratio (T n (t)/T n -i(t)) [ll[ or sufficiently large energy 
separation A/x, in which case r„(t^)/r„_i(i*_ 1 ) can be 
large even if T T1 (t) w r„_i(i) due to the difference in t b n _ l 
and t b n . Alternatively, devices with reduced coupling be- 
tween barrier and plunger may benefit from the thermal 
limit which at sufficiently low temperature will further 
enhance the initialization precision. This regime may be 
reached by adding compensation pulses to Vgd during 
the transition (i) — > (ii) in Fig. [2] This example illus- 
trates the importance of identifying the capturing mech- 
anism for an efficient optimization of the initialization 
process of dynamic QDs as building blocks for nanoelec- 
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Appendix: supplemental material on data evaluation 

Here we detail the extraction of probabilities for the 
number of electrons n which are captured on the quantum 
dot, subsequently transferred onto the node and observed 
by the detectors labeled Detl and Det2 in Fig. 1 in the 
main text. 

In contrast to other detection methods, like e.g. quan- 
tum point contacts, the superconducting SETs used here 
to monitor the node's charge state respond periodically 
to external potentials thus showing a specific response 
(amplitude as well as direction) depending on the SET's 
working point. This is influenced by all potentials in close 
vicinity to the SET island via capacitive coupling. In ad- 
dition to the gates controlling the semiconducting part of 
the chip by applying voltages Vqs , Vgd and V" barrior we 
have evaporated gates close to each of the SET islands 
(referred to in the following as SET gates with voltages 
Vgset) which allow for the control of the SET working 
points independently from any other potential in the cir- 
cuit. The gate voltage Vgset is always the same for both 
SET gates. A typical current response for both detectors 
under variation of Vgset is shown in Fig. IS5I For each 
value of V"gd (gate between the quantum dot and the 
node), we change the voltage Vgset (and thereby the 
SET working points) in 20 steps over half a SET period. 
The parameter range is indicated by the black vertical 
lines in Fig. IS5l which corresponds to Vgset € [—2.5,5] 
mV. 

For each setting of the static gates V*gd and Vgset we 
now take a timetrace of the detectors' current response 
to a specific pulse sequence (see Fig. [S6]) : We first ap- 
ply a short pulse to the gate Vbarricr between the node 
and the drain lead (black vertical line) which randomly 
changes the number of electrons stored on the node close 
to an equilibrium number. Due to this pulse, we avoid 
the buildup of a strong electrostatic potential on the node 
which on the one hand would couple capacitively to the 
quantum dot (mesoscopic feedback) but would also po- 
tentially drive the SETs out of sensitivity. After this 
first pulse, three combined capture and transfer pulses 
onto gate Vgs of the dynamic quantum dot are applied 
which arc illustrated by green vertical lines in Fig. IS6I 
All pulses are delayed by fdeiay = 40 ms each. Every 
timetrace consists of 279 sequence periods, i.e. 279 reset 
pulses and 837 combined capture and transfer pulses. 

Next we describe an example of evaluating the trace 
shown in Fig. IS6I We first make a histogram (see Fig. [S7l 
and lS8p of the first 175 sequence periods of the two time 
traces, then apply a second-order Savitzky-Golay filter 



in order to smooth the data and calculate the derivative 
(black lines), which yields the maxima and minima of the 
histograms of the time traces. The maxima found cor- 
respond to the occupation levels of the node. The data 
range of the SETs indicating the range of sensitivity is 
also derived from this calculation and is limited by the 
first and last minimum found by the algorithm. These 
either mark the vanishing separation of levels close to 
the SET extrema (lower limit for Detl) or possibly in- 
sufficient occupation of the following levels for reliable 
detection (upper limit for Detl, both limits for Det2). 

The resulting levels are then compared to each mean 
value of the corresponding detector trace between two 
pulses. Based on the assumption that each possible 
charge state of the node is found by the level detec- 
tion (which is supported by the stochastic reset of the 
node's charge number in each sequence) we can directly 
extract the number of charge carriers between two pulses 
by counting the number of levels lying in-between the 
starting and the final level. Whenever the starting or 
the final state is beyond these limits set by the algorithm 
(i.e., the SET is not sensitive) the corresponding event is 
not considered in the statistics. 

This evaluation is performed for every parameter set 
(Vgd, Vgset)- The results are combined for such traces 
where Vgd is constant, both of the detectors are sensi- 
tive for at least 30% of the time and more than 150 equal 
events for both detectors have been observed. The last 
argument is based on the fact that both detectors are 
coupled to the same island and therefore the two signals 
are not independent but correlated in terms of electron 
number difference between two pulses. By applying this 
restriction we can strongly reduce the number of erro- 
neous counts caused by noise on one of the detectors or 
by failures in the algorithm (mainly errors in the level de- 
tection). All events or traces not matching these criteria, 
e.g. pulses at insensitive working points for at least one 
SET or inequal results for both detectors, are neglected. 
Therefore, the total number of events A(Vgd) evaluated 
for one setting of Vgd varies from 1853 to 8997 events 
for Vgd 6 [—230,-175] mV (average number of events 
« 6100). 

Finally, this procedure results in a sequence of electron 
numbers captured and transferred by the pulses applied 
to the dynamic QD. For each Vgd we count the num- 
ber of outcomes M„(Vgd) that have resulted in exactly 
ri = ... 4 electrons captured and transferred in one cy- 
cle. M„(V"gd) is expected to obey the binomial distri- 
bution with the (unknown) success probability P n (VQu) 
and the total number of trials -/V(Vgd)- The probability 
P n is then estimated by definition as success frequency: 
P n = M n /N. We have used the principle of maximal like- 
lihood and the quantile function of the binomial distri- 
bution to estimate the (symmetric) p = 95% confidence 
intervals for P n , as depicted by error bars in Fig. 3 and 4. 
For the special outcomes of M = and M = N the con- 
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fidence interval is estimated by requesting the expected 
probability of this extremal outcome to be no less than 
P- 

For the best-fit estimation of parameters of the mod- 
els expressed by Eq. (6) and (9a) or Eqs. (8) and (9b) 
we have used a standard least-squares routine, matching 
the Gaussian model it assumes to p — 68% confidence 
intervals obtained from the exact binomial statistics as 
described above. The confidence intervals for the model 
parameters are then computed using Monte Carlo simula- 
tion of 500 synthetic data sets for each of the two models 
separately. The synthetic data are drawn from binomial 
distributions with success probabilities P„(Vgd) taken 
from the respective best-fit model and the number of tri- 
als .ZV(Vgd) equal to that in the actual experiment for 
each Vqd- 




Figure S5: Typical SET responses under variation of the SET 
gates. The red trace corresponds to Detl, the blue one to 
Det2. The vertical black lines indicate the range of gate volt- 
age settings used for working point variation of the SETs. 




Figure S7: Histograms for both traces together with the first 
smoothed derivative and the maxima/minima found indicated 
by the dots. The vertical dashed lines show the limits of 
sensitivity set by the first and last minimum found. 
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Figure S8: Signal traces together with the histogram and the 
levels found (dashed horizontal lines in the trace) and the 
range of SETs' sensitivity (bold dashed lines). 
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Figure S6: Signal traces for Detl (lower graph, red trace) 
and Det2 (upper graph, blue trace) both measuring the same 
signal. Vertical lines indicate pulse positions as explained in 
the text. 
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